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Abstract 

Experiments involving copper-copper collisions at the RHIC-BNL (USA) at energies 
y/sN~N = 22.5, 62 and 200 GeV have produced a vast amount of high-precision data which 
are to be analysed in the light of various competing models in the domain of multiparicle 
production scenario. We have chosen to analyse here the measured data on the py-spectra of 
various light and non-strange secondaries at various energies mentioned above, some of their 
very important ratio-behaviours at the various centralities of the collisions and the nuclear 
modification factors, Raa and Rcp, in the light of a version of the Sequential Chain Model 
(SCM). The agreements between the measured data and model-based results are generally 
found to be modestly satisfactory. Besides, our obtained results have also been compared 
with results based upon two approaches with strong standard model flavour, of which one is 
purely pQCD-oriented and is forwarded by Vitev. 
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1 Introduction 



In a previous paper [T] we dealt with the properties of Au + Au interactions at RHIC-energies at 
y / sjv r /v=200 GeV in great detail with emphasis on the nature of the most important observables. 
The gold-gold interactions constituted a relatively heavy system. Compared to it, the Cu + Cu 
system represents a lighter and smaller system-size, with mass-number for copper being much 
less than gold. Besides, in Ref. pQ we also gave the details about the essential points of the 
basic multiparticle production model applied therein. 

So, in order to avoid repetitions, we will just give a very brief sketch of the main/major 
observables that would here be dealt with and on the same model that we make use of here. 

Our purpose is to check (i) whether the same model could explain the data measured for a 
comparatively lighter system; (ii) whether [and to what extent (if at all)] there is any prominent 
system size effect on the results obtained by the experimental measurements. 

Apart from the nature of p^-spectra and some particle-production ratios, we will dwell upon 
here properties of nuclear modification factors Raa and Rcp for Cu + Cu collisions and the 
centrality dependence of the results. But production of the heavies and the heavy strange 
hadrons are left out in the present work. 

As our approach to the analysis of Cu + Cu collisions (or for that matter, of Au + Au 
and even p + p reactions) is somewhat different from the 'standard' varieties, our phrases and 
vocabularies are not identical with them. We have inducted here the physics of partonic multiple 
scatterings for both nucleus- nucleus interactions, and p + p reactions at large-pr through some 
mechanism as given in the text in the appropriate place. Even the effects of rescatterings have 
also been incorporated here by the same technique. But, the associated hadronic collectivity 
factor (gauged by the elliptic flow represented normally by 'iV) has not been touched upon here; 
this flow-behaviour would be dwelt-upon in some detail in a future work. 

In order to dispel any misunderstanding of any interested reader, let us emphasise here, before 
closing up the introductory remarks, that our intention is not to make any tall claim that the 
model tried out and tested here is much better than the existing models. But, our objective 
is just to familiarise the fact that modest description of data is quite possible even with some 
non-standard approaches, like the present one, which is being applied and validated to explain 
some observables in high energy collisions since the mid-seventies. 

Finally, the organization of this work is as follows. In section 2 we sum up the main physics- 
aspects entailed in this work. A glance and a glimpse into it would reveal the non-standard 
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nature of the basic approach. In the next section (section 3) the results arrived at have been 
presented with tables and figures. The section 4 provides some comparison of our work with 
both data and the results obtained by other model-based approaches for a very few observables 
related to the production of some very select secondaries (for which theory-based studies on 
Cu + Cu collisions are available). And in the last section (section 5) we offer the final comments 
and conclusions. 

2 The Theoretical Framework: A Brief Outlook 

The description of the model-based features would be subdivided into some parts. The first 
part (subsection 2.1) gives a brief overview of the production mechanism of the secondary hadrons 
in nucleon-nucleon (p+p) interaction in the context of the Sequential Chain Model (SCM). Then 
in the subsection 2.2 we present a brief outline of the main and major achievements of the model; 
these points, in essence, also highlight the important characteristics of the model. Thereafter, 
the relevant transitions for different observables from p+p to A + B interactions will be discussed 
in the subsection 2.3. 

2.1 Basic Model for Particle Production in PP Scatterings: An Outline 

According to this Sequential Chain Model (SCM), high energy hadronic interactions boil down, 
essentially, to the pion-pion interactions; as the protons and neutrons are conceived in this 
model as p = (7r + 7r i9) and n = (1/ \Jx\ + x| [Xi( 7r ° 7r °'#)+X2(7r ~7r + $)]) repectively, where ■& is 
a spectator particle needed for the dynamical generation of quantum numbers of the nucleons 
and xii Xi are t ne weightage factors PQ- [5]. Our focus and concern would now be confined 
and concentrated to the structure of protons alone. The production of pions in the present 
scheme occurs as follows: the incident energetic ir- mesons in the structure of the projectile 
proton (nucleon) emits a rho(^)-meson in the interacting field of the pion lying in the structure 
of the target proton, the ^-meson then emits a 7r-meson and is changed into an omega(w)-meson, 
the u;-meson then again emits a 7r-meson and is transformed once again into a ^-meson and thus 
the process of production of pion-secondaries continue in the sequential chain of q-co-tt mesons. 
The two ends of the diagram contain the baryons exclusively [I]-[5]- 

For K + {K~)ox K°K° production the model proposes the following mechanism. One of the 
interacting 7r-mesons emits a ^-mesons; the g-mesons in its turn emits a ^°-meson and a 7r-meson. 
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The 7r-meson so produced then again emits g and <fP mesons and the process continues. The 
mesons so produced now decays into either K + K~ or K°K° pairs. The q-tt chain proceeds in 
any Fenymann diagram in a line with alternate positions, pushing the cfP mesons (as producers 
of K + K~ or K°K° pairs) on the sides. This may appear paradoxical as the (fP production cross- 
section is generally smaller than the KK production cross-section; still the situation arises due 
to the fact that the cfP resonances produced in the collision processes will quickly decay into KK 
pairs, for which the number of <jp will be lower than that of the KK pairs. Besides, as long as </>° 
mesons remain in the virtual state, theoretically there is no problem, for <fi°K + K~ ( or (j) K°K ) 
is an observed and allowed decay mode, wherein the strangeness conservation is maintained 
with the strange-antistrange coupled production. Moreover, (fpK + K~ ( or (p K°K ) coupling 
constant is well known and is measured by experiments with a modest degree of reliability. And 
we have made use of this measured coupling strength for our calculational purposes, whenever 
necessary. It is assumed that the K + K~ and K°K° pairs are produced in equal proportions 
[1J-|5J. The entire production process of kaon-antikaons is controlled jointly by the coupling 
constants, involving g-ir-cft and (fi°-K + K~ or (j) -K°K . 

Now we describe here the baryon-antibaryon production. According to the SCM mechanism, 
the decay of the pion secondaries produces baryon-antibaryon pairs in a sequential chain as 
before. The pions producing baryons-antibaryons pairs are obviously turned into the virtual 
states. And the proton-antiproton pairs are just a part of these secondary baryon-antibaryon 
pairs. In the case of baryon-antibaryon pairs it is postulated that protons-antiprotons and 
neutrons-antineutrons constitute the major bulk, Production of the strange baryons-antibaryons 
are far less due to the much smaller values of the coupling constants and due to their being much 
heavier. 

The field theoretical calculations for the average multiplicities of the tt, K and p-secondaries 
and for the inclusive cross-sections of those secondary particles deliver some expressions which 
we would pick up from pQ-[5]. 

The inclusive cross-section of the 7r _ -meson produced in the p + p collisions given by 

E^t\ pp ^- x ^ IV exp(-2.38 < n„.- > pp x) 1 exp(- ~ A ^ I>J ■ < 1 < 



with 



dp 6 pp (N% ) < > pp (1 

Pt 



< n n + > pp ^ < TV- >p P = < n n o > pp l.ls 1/5 , (2) 



where r_- is the normalisation factor which will increase as the inelastic cross-section increases 



and it is different for different energy region and for various collisions, for example, \T v -\ = 90 
for Intersecting Storage Ring(ISR) energy region. The terms pr, x in equation (1) represent 
the transverse momentum, Feynman Scaling variable respectively. Moreover, by definition, 
x = 2pi/ yfs where pl is the longitudinal momentum of the particle. The s in equation (2) is 

the square of the cm. energy. 

N w ~ 

l/p T R of the expression (1) is the 'constituent rearrangement term' arising out of the partons 
inside the proton which essentially provides a damping term in terms of a power-law in px with 
an exponent of varying values depending on both the collision process and the specific pr-range. 
The choice of Nr would depend on the following factors: (i) the specificities of the interacting 
projectile and target, (ii) the particularities of the secondaries emitted from a specific hadronic 
or nuclear interaction and (iii) the magnitudes of the momentum transfers and of a phase 
factor (with a maximum value of unity) in the rearrangement process in any collision. And 
this is a factor for which we shall have to parameterize alongwith some physics-based points 
indicated earlier. The parametrization is to be done for two physical points, viz., the amount of 
momentum transfer and the contributions from a phase factor arising out of the rearrangement 
of the constituent partons. Collecting and combining all these, we proposed the relation to be 
given by [6] 

N R = 4 < N part > x / 3 8, (3) 

where < N part > denotes the average number of participating nucleons and 9 values are to 
be obtained phenomenologically from the fits to the data-points. In this context, the only 
additional physical information obtained from the observations made here is: with increase in 
the peripherality of the collisions the values of 6 gradually grow less and less, and vise versa. 
Similarly, for kaons of any specific variety ( K + , K~, K° or K° ) we have 

d?a 1 —1 33p 2 

E-ra\pp-+K-<B ~ T K- exp(-6.55 < n K - > pp x) exp( : , (4) 

a P J n r > <Uu- > pp 



Pt 



with |IV-| = 11.22 for ISR energies and with 



< n K + > pp =< n K - > pp =< n K o > pp =< n R(s > pp ^ 5 x 10 V /4 . (5) 

And for the antiproton production in pp scattering at high energies, the derived expression for 
inclusive cross-section is 



5 




(6) 



< Up > 



=< n v > 



pp 



2 x 1(T 2 s 1 ^ , 



(7) 



2.2 Some Cardinal Characteristics and Triumphs of the Model 

We agree that the word 'non-standard' used in the preceding section is a port manteau adjective, 
with many layers of meaning hidden within it. In this particular case, (i) the model is based on 
some new ideas about the structure of hadrons and the nature of hadronic interactions; (ii) the 
proposed mechanism underlying this work does not admit of any airtight compartmentalisation 
of the 'soft' (low-pr) and 'hard' (large pr, Pt > 2 GeV/c) production; (iii) rather, the model 
presents a unified approach to the production of particle-secondaries; (iv) besides, the fundamen- 
tal expressions for final (analytical) calculations are derived here on the basis of field-theoretic 
considerations and the use of Feynman diagram techniques (albeit with some simplifying high 
energy approximation and assumptions) with the infinite momentum frame tools and under 
impulse approximation method; (v) this approach establishes the < universality > aspect of 
the multiplicity of high energy interaction called 'globality property'; (vi) the model explains 
the 'jet '-structure for emanation of the secondary particle as the << two — sidedsprays >> of 
hadrons, (vii) it reproduces the behaviour of average multiplicity, nature of (invariant) inclusive 
cross-section and the properties of average transverse momenta of various secondaries; (viii) the 
model could also account for the very slow rising nature of the total cross-sections, (ix) Besides, 
this model can/does explain the majority of the characteristics of what are known or believed 
to be the 'quark gluon plasma'-diagnostics. And that could be obviously done by an alternative 
approach and outlook, (x) Furthermore, the application of the model can also accommodate a 
large amount of very important Cosmic Ray Physics issues and problems that came to the fore 
in the very recent times. 

The above-mentioned features are essentially ingrained in the entirety of this totally non- 
standard mechanism emerging from an alternative philosophy and outlook about the particle- 
constituents and their interaction mode. On the whole, this is purely an analytical approach with 
a reasonable number of valid assumptions and approximations which are commonly used by all 
High Energy Physicists. Uptil now, we have confined ourselves mostly to the non-simulational 
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calculations. The calculations for 'soft'and 'hard' have been superposed here by virtue of the 
simple factorisation property. One of the very strong points about this model is the fact that the 
various coupling strengths used in this model are not only known but also reliably well-measured 
by several experiments. This factors helps to reduce considerably the speculative components 
in the results. 

2.3 Results for AA Collisions from PP Reactions: The Connecting Bridge 

In order to study a nuclear interaction of the type A + B — >• C~ + x, where A and B are projectile 
and target nucleus respectively, and C~ is the detected particle which, in the present case, would 
be ir~ , K~ and p, the SCM has been adapted, on the basis of the suggested Wong [7] work to 
the Glauber techniques by using Wood-Saxon distributions [8j-[9]. The details of calculations 
and the features of the SCM have been given in our previous paper pQ. 

The general form of our SCM-based transverse-momentum distributions for A + B — > C~ + X- 
type reactions can be written in the following notation: 

a+b-,c-+x = ol c - exp(-/3 c - x p T ). (8) 



27TPT dprdy^^ ^ ^ no 

The set of relations to be used for estimating the parameter olq- is given below [lj. 

(Aa B + Ba A ) 1 

a °~ = ^ l + a{A^ + B^) Tc - eM ~ V < nc ~ >pp X) (9) 

Here, in the above equation [eqn.(9)], r^-, as stated above, is the normalization constant which 
is different for the different secondaries and the collider energies. It also depends on the centrality 
of the collisions. The first factor in eqn.(9) gives a measure of the number of wounded nucleons 
i.e. of the probable number of participants, wherein Agb gives the probability cross-section of 
collision with t B > nucleus (target), had all the nucleons of A suffered collisions with i?-target. 
And BaA has just the same physical meaning, with A and B replaced. Furthermore, a a is 
the nucleon(proton)-nucleus(A) interaction cross-section, <jb is the inelastic nucleon(proton)- 
nucleus(B) reaction cross-section and <jab is the inelastic AB cross-section for the collision of 
nucleus A and nucleus B. The values of gab-, &a, &b are worked here out in a somewhat heuristic 
manner by the following formula [TO] 



<ffi = ^ (AH 3 ojecUle + Al£ get - 5f (10) 



with cr = 68.8 mb, 5 = 1.32. 



Besides, in expression (9), the second term is a physical factor related with energy degradation 
of the secondaries due to multiple collision effects. The parameter a occurring in eqn.(9) above 
is a measure of the fraction of the nucleons that suffer energy loss. The maximum value of a 
is unity, while all the nucleons suffer energy loss. This a parameter is usually to be chosen [7], 
depending on the centrality of the collisions and the nature of the secondaries. 

The values of r\ in eqn. (9) are different for different secondary produced; for example, rj= 
2.38 for pions, 6.35 for kaons and 25.4 for protons, as were given in eqn.(l), eqn. (4) and eqn. (6). 

N c ~ 

l/p T R of the expression (8) is the 'constituent rearrangement term' arising out of the partons 
inside the proton which essentially provides a damping term in terms of a power-law in px with 
an exponent of varying values depending on both the collision process and the specific p^-range. 
We have already mentioned the details earlier. 

The values of (3q- of the equation (8) for different secondaries have been calculated with the 
help of eqn.(l), eqn. (2), eqn. (4) -eqn. (7). 

3 Steps Towards Calculations 

At the very start let us present a summary (which might be a rehash of what we have men- 
tioned) of the key physical facts that would be of paramount importance which proceeding 
towards calculations. There are some foundational steps that enable us to arrive at the final 
working formulae which deliver the results to be reported here for for Cu + Cu collisions. The 
procedural steps are as follows: (i) Firstly, we have the basic model for p + p scattering at high 
energies and low-pT ('soft') interactions, so we have to convert the mathematical expressions 
for nucleus+nuclus (Cu + Cu) collisions by introducing the nuclear dependence factor on the 
results arrived at for p + p collisions, (ii) Secondly, the data-points on Cu + Cu reaction at 
various high energies exceed the range of the low-p? boundary, pr > 2 GeV/c, for which the 
large-pT effect is to be superposed on the expressions for soft-production of the secondaries. The 
constituent (partonic) rearrangement factor introduced here takes care of this physical feature 
arising out of the 'hard' (large-pr) contributions, (iii) In a model-dependent way the SCM has 
some special and specific means of excess production of the positive secondaries, specifically the 
light secondaries. 
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Table 1: Values of (a 7r -) pp , (N^ ) pp and (f3 w -) pp for it productions in p + p collisions at 
^s7w=20, 63 and 200 GeV 





( a iT- )pp 


(Nf)pp 


(Ar~ )pp 


20 GeV 


0.281 


4.086 


0.703 


63 GeV 


0.545 


3.327 


0.468 


200 GeV 


0.007 


3.867 


0.293 



3.1 Production of Main Varieties of Negatively Charged Secondaries 

The general expressions of inclusive cross-sections for the production of ir~ , K~ and p for p + p 
collisions were stated in the previous section by eqn.(l), eqn.(4) and eqn.(6) respectively. 

For the production of 7r~-mesons m.p + p collisions, we use eqn.(l), eqn.(2) and eqn.(9) with 
a n - = T n - exp(— 2.38 < n n ~ > pp x). The values of (a n -) pp , (iVjj ) pp and {(3 7T -) PP are given 
in Table 1. The experimental data for the inclusive cross-sections versus pT[GeV/c] for ir° 
production in p + p interactions at y / J/v r /v = 20 GeV are taken from Ref. [12] . And for the data 
for inclusive cross-sections for ir + at energies ^snn = 63 GeV and ^/snn = 200 GeV we use 
references [13], [13] , (1~5] respectively. They are plotted in Figure 1(a), Figure 1(b) and Figure 
1(c) respectively. The solid lines in those figures depict the SCM-based plots. 

3.1.1 Production of 7r -mesons in Cu + Cu Collisions 

We now, at first, turn our attention to tt~ production in Cu + Cu collisions at energies ^snn 
=22.5, 62.4 and 200 GeV. 

Using eqn.(8) and eqn.(9) the SCM-based expressions for transverse momentum distribution 
of negative pions produced in the Cu+Cu collisions at y / i]vjv = 22.5, 62.4 and 200 GeV at RHIC 
and for different centralities can be obtained. The values of a w - , and f3 n - for different 

centralities and for different energies are given in Table 2. The values of N par t, for calculating 
Nft from eqn.(3), in this context, have been taken from Ref. [16], |17j . The experimental 
results for the production of 7r~ at different centralities for energies y / S7v r /v = 22.5, 62.4 GeV are 
taken from the Ref. [18] and for energy ^snn = 200 GeV we have used Ref. [19] . The invariant 
yields for tt~ against pT[GeV/c] for different energies are plotted in Figures 2(a), 2(b) and 2(c) 
respectively. The solid lines in those figures show the theoretical SCM results. 
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3.1.2 K Production in Cu + Cu Collisions 

With the help of the eqn.(3)-eqn. (5), eqn. (8) and eqn. (9) the values of a^-, Nj£ and Pk- 
for the transverse momentum distributions for different centralities of ^"-particles in Cu + Cu 
collisions at -y/s/viv = 22.5, 62.4 and 200 GeV at RHIC have been calculated and they are given 
in Table 3. The experimental results are taken from the PHENIX group [18] for RHIC energies 
^snn = 22.5 and 62.4 GeV and they are plotted in Figures 4(a) and 4(b), whereas for y / J/vjv 
= 200 GeV, the Ref. [19] has been used and data are plotted in Fig. 4(c). The lines in those 
figures depict the theoretical outcomes. 

3.1.3 Production of p in Cu + Cu Collisions at Different Energies 

Using the eqn.(6)-eqn.(9), the transverse momentum distributions for antiproton in Cu + Cu 
collisions at energies ^snn = 22.5, 62.4 and 200 GeV at RHIC have been calculated the cor- 
responding values of ap, and (3p for different centralities and different energies are given in 
Table 4. The experimental results of invariant yields for the production of antiproton produc- 
tions for energies ^Jsnn = 22.5, 62.4 GeV and for ^Jsnn = 200 Gev are taken from the Refs. 
[18] . [19] respectively. They are plotted against px in Figures 6(a), 6(b) and 6(c) for different 
centralities. The solid lines in those figures show the theoretical SCM-based results. 

3.2 On Excess Production of Positive Particles 

True, on the average, the particles are produced in charge-independent equal measure, for which 
roughly one-third of the particles could be reckoned to be positively charged, one third are 
negatively charged and the rest one third are neutral. But, according to the present mechanism 
of particle production, there are some specifically exclusive means to produce positive particles of 
which 7r + , if , and p are the members. They are produced from within the structure of protons 
(nucleons). These production characteristics and the quantitative expressions for their special 
production have been dwelt upon in detail in Ref. [5]. Let us assort the relevant expressions 
therefrom as results to be used here. 

3.2.1 Production of Positive Pions in Cu + Cu Collisions 

For production of positive pions [ir + mesons] the excess term could be laid down by the following 
expressions [5]: 
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where the symbols have their contextual connotation with the following hints to the physical 
reality of extraneous ir + , as non-leading secondaries. The first parts of the above equations 
(Eqn.(ll)), contain the coupling strength parameters, the second terms of the above equations 
are just the propagator for excited nucleons. The third terms represent the common multiparticle 
production amplitudes along with extraneous production modes and the last terms indicate 
simply the phase space integration terms on the probability of generation of a single tt + . These 
expressions are to be calculated by the typical field-theoretical techniques and are to be expressed 
- if and when necessary - in terms of the relevant variable and/or measured observables. 

In order to arrive at the transverse momentum distribution of 7r + , one has to consider the 
Eqn. (1), eqn. (8) along with eqn. (11). For excess tt + production, a factor represented 
by (1 + 7 7r+ p^ + ) is to be operated on ^:£^\cu+Cu^-+x as an multiplier 0. i* + ~ 
(207r$^/ < n-x >)/y/s ~ 0.44 PQ, [IE]. Taking < p T > 7r +^ 0.31 GeV/c [20], the calculated 
values of a n + , and /3 7r + for different centralities and for different energies are given in Table 
2. In Figures 3(a), 3(b) and 3(c), we have plotted experimental versus theoretical results for 
7r + production in Cu + Cu collisions at energies ^snn = 22.5, 62.4 and 200 GeV, respectively. 
Data are taken from the Refs. [18] and |19| . The solid lines in those Figures are the SCM-based 
plots. 

3.2.2 Production of Positive K-mesons in Cu + Cu Collisions 

For the excess production of K + -mesons we proceed in the same path as we did earlier for the 
case of 7r + production. The equation for the extraneous production of K + is given hereunder 

1 If d 3 kx 

(B k +)pp = ^9kna + ^91kn) ~ R y _ m 2}2 A ( u > J 2fc (27r) 3 ex PH fc ^ x )' ( 12 ) 

Adopting the above procedure, as we indicated for the production of positive pions, we obtain 
for the transverse momentum distribution of K + a multiplicative factor ~ (1 + 7^ + p^ + )to be 
operated on d%>Tdy I Cu+Cu-^K- +x 0- For the production of K + , the factor, calculated from 
eqn. (12), ~ {^9knk + ^ 7r g1 KN )/2^/s ~ 0.082 [I], [TT]- We use the value of < p T >k+- 
0.36 GeV/c [20] and the corresponding values of ax+, + and (3 K + for different energies are 
presented in the Table 3. The experimental results [18], [19] for the production of K + of different 
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centralities and for different energies are plotted in Figures 5(a), 5(b) and 5(c). The solid lines 
in those figures show the theoretical plots. 

3.2.3 Production of Excess Protons 

Similarly, for the excess production of protons, the extraneous term can be can be picked up 
from our previous work [5] in the following form: 



For the production of protons, we obtain for the transverse momentum distribution of p by 
operating a multiplicative factor ~ (1 + 7 p p^), which is an outcome of eqn. (13), on 2 ^ PT dp^dy • 
The value of 7 P ~ 0.32 [1J and by taking < p T > p ~ 0.50 GeV/c [20], we finally obtain the 
values of a p , and (3 P , which are given in the Table 4. In Figures 7(a), 7(b) and 7(c), we have 
presented the experimental values of invariant yields for proton-production versus the theoretical 
SCM-based results for energies y / J/vjv = 22.5, 62.4 and 200 GeV respectively. Data are taken 
from [18] and [19]. The lines in the figures show the theoretical outcomes. 

3.3 The Ratio Behaviours for Different Secondaries 

3.3.1 The vr /vr+ Ratios at ^/s^ = 62.4 and 200 GeV 

The model-based 7r~ /ir + ratios for different participating nucleons N part at energies s/snn = 
62.4 and 200 GeV have been obtained from the expression (8) and Table 2. Data in Figs. 8(a) 
and 8(b), shown by filled squares and blank circles respectively, are taken from the PHOBOS 
group [21], [22]. The theoretical values in this regard are plotted by solid line in Figure 8(a) and 
by filled circles in 8(b). 

3.3.2 The K~/K + Ratios at ^/sW^ = 62.4 and 200 GeV 

In a similar way, the N par t versus K~ /K + at y / J/vjv = 62.4 and 200 GeV can be obtained from 
equation (8) and Table 3. The calculated values of K~ /K + against the N par t in the light of 
SCM are shown by solid line in Fig. 9(a) at energy ^snn = 62.4 GeV and by solid squares in 
Fig. 9(b) at ^/sWn = 200 GeV. The data in those figures are taken from PHOBOS |2T]. [22]. 




(13) 
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3.3.3 Some Other Ratio-Behaviours at ^/s^ = 62.4 and 200 GeV 

Based on the SCM, the p/p ratios at different energies like ^/snn = 62.4 and 200 GeV and for 
different participating nucleons N par t are obtained with the help of equation (8) and Table 4. 
The calculated values are plotted in Figs. 10(a) and 10(b) at energies 1 /5aF/v = 62.4 and 200 
GeV by solid line and filled circles respectively. Data in those figures are taken from PHOBOS 
EH, [22]. 



3.3.4 The p/vr~ and p/ir + Ratios at ^/sWn = 22.5, 62.4 and 200 GeV 

The expressions for p/ir~ ratios against pt for central reactions at energies ^snn = 22.5, 62.4 
and 200 GeV can be obtained from equation (1) and Table 2 and Table 4 and they are given 
hereunder 

— = 0.32p^ 038 exp(-0.15^) for V5^7 = 22.5 GeV, (14) 

7T 

■4: = 0.28^ 932 exp(-0.15^) for ^^ = 62.4 GeV, (15) 

7T 

■4- = 0.22p^ 821 exp(-0.13p|) for ^1^ = 200 GeV. (16) 

And for p/ir + ratios versus pr at energies ^snn = 22.5, 62.4 and 200 GeV the SCM-based 
equations are written as 

-4 = 0.85^ 038 exp(-0.15?4) for ^1^ = 22.5 GeV, (17) 

-4 = 0.42^ 932 exp(-0.15p^) for = 62.4 GeV, (18) 

-4 = 0.35p^ 821 exp(-0.13p^) for = 200 GeV. (19) 

In Figures 11(a) and 11(b) we have plotted px versus p/vr _ and p/n + respectively for central 
Cu + Cu collisions at ^Jsnn = 22.5, 62.4 and 200 GeV. Data of these figures are taken from 
PHENIX [23J. Lines in Figures 11(a) and 11(b) represent eqn. (14)- eqn. (16) and eqn. (17) to 
eqn. (19) respectively. 

3.4 Nuclear Modification Factors 

In this subsection we would dwell upon the Nuclear Modification Factors of two types viz., Raa 
and Rcp- In 3.4.1 the former (Raa) would be defined and the results would be hinted, though 
the figures for both would be shown in the next section. And in 3.4.2 the second one (Rcp) 
would be treated in some detail. 
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3.4.1 The Nuclear Modification Factor, Raa 

The nuclear modification factor (NMF), designated as Raa-, for any secondary C, is denned by 

m 

R c {l/N%j)<PNC A /dp T dy 
AA < Ncottf) > x tPog/dprdy- { ' 

Depending on this definition, the SCM-based results on NMFs for Cu + Cu collisions at y / J/vjv 

= 22.5, 62.4 and 200 GeV are deduced on the basis of Eqn.(8), Table 1 and Table 2 and they 

are given by the undernoted relations 

R AA = 0.730p^ 621 for ^/J^ = 22.5 GeV, (21) 

Raa = 0.680^ 311 for ^/J^ = 62.4 GeV, (22) 

Raa = 0.315p£ 234 for y/wfi = 200 GeV. (23) 

wherein the values of < N co u(b) > to be used are « 140.7, 152.3 and 182.7 [17J for Cu + Cu 
collisions at y/s^N = 22.5, 62.4 and 200 GeV respectively. For N evt , we use the values « 
5.8 x 10 6 , 192 x 10 6 and 794 x 10 6 [T7] at three different energies and a pp to be used as 30 mb 
. Our model-based plot is shown in a figure in the next section. 



3.4.2 The Nuclear Modification Factor, R C p 

There is yet another nuclear modification factor, Rcp which reflects precisely the hadron 
spectra in different centrality bins and presents comparison of the py-spectra between a collision 
at a specific centrality and the one at the relatively peripheral collision. It is quantified as 



[d 2 N/(2np T dp T dy)/N bi n 



central 



R ° P [d 2 N/(2irp T dp T dy)/N bm ]P en P heral ' ^ 
According to above definition, with the centrality set at 0—10% and the peripherality at 60—94%, 
our model-based expression for Rcp in case of neutral pions (7r°s) is given by 

Rgp = 0.543^ - 054 . (25) 

wherein we have made use of the values shown in Table-2. The Rcp for neutral kaons (K°s) 
with the same centrality and the peripherality changed to 40 — 60%, our model-based result is 

R%°p = 0.664^ a235 . (26) 

wherein we have taken K° = 1/2(K + + K~), and used the necessary values from Table-3. 
Our model-dependent plots on these two particular observables, Raa and Rcp-, are shown in 
the subsequent section on comparative studies. 
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4 Data and Results on Some Select Observables: A Comparison 
between Models 

We strongly uphold the view that point-to-point or secondary-to-secondary specific comparisons 
between our model-based results with both data and the other model-based calculations would 
be quite meaningful and physically significant. But the main and major constraint in this 
attempt is the lack of availability of such comprehensive calculations encompassing all the light 
secondaries. In so far as Cu+Cu collisions are concerned, we have, so far, come across two model- 
based studies, of which one is the application of the Quark Combination Model (QCM) made 
by Fei et al. [27] for the production of neutral pion/kaon and the other is the pQCD-oriented 
theoretical study done by Vitev |28j for production of only the neutral pion and for no other 
secondary. This constitutes a gross limitation to the successful completion of the comparison- 
aspects in the present study. However, we have tried here to show some comparison(s) only for 
one or two observables related to a few selected neutral secondaries, with whatever little other 
model-dependent studies could be obtained uptil now. 

In spite of the difficulties mentioned in the above paragraph, in the adjoining Fig. 12 we 
compare the data-versus-results based on two models for production of (a) neutral pion and (b) 
neutral kaon m. Cu + Cu collision at y / ijv r /v = 200 GeV. The standard variety of models that 
are reckoned with here for comparison with the present SCM are (i) The Quark Combination 
Model (QCM) and (ii) Perturbative QCD-inspired Vitev's Model. In the former (QCM), the 
main idea is to line up the quarks and antiquarks in a one-dimensional order in phase space, 
e.g. in rapidity and let them combine into initial hadrons one by one following a combination 
rule [27] . These initial hadrons through combination of constituent quarks are then allowed to 
decay into the final state hadrons through the decay program of PYTHIA 6.1. The calculational 
steps, in this model, proceed on the basis of two-prong assumptions of two-component model 
(based on 'soft'-'hard' artifact) and the concept of parton-hadron duality. 

The latter model used by us for comparison is one of Vitev [28]. This is essentially a pQCD- 
oriented model with an analytic model of jet-quenching which embraces medium-induced energy 
loss after hard partonic scattering. This approach reduces the jet cross-section in the presence 
of the medium but leaves the parton fragmentation function unaltered. In fact, this feature 
was conveniently implemented in the analytic model of Vitev [28] for QGP-induced leading 
hadron suppression. Vitev actually made use of a Hagedornian form of power law expressions 
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for invariant cross-section(s) and inducted also the radiative energy-loss formalism. 

In Fig. (12), the invariant yields versus px (GeV/c) for (a) tt° and (b) for different 
centralities in Cu + Cu collision at ^snn = 200 GeV have been plotted. Data are taken from 
[29j and [30]. The solid lines in those Figures represent the SCM-based calculations wherein the 
dashed lines show the QCM-oriented results [27]. And the dotted line in Fig. 12(a) shows pQCD- 
inspired calculations done by Vitev [28] for production of only the neutral pions in Cu + Cu 
collision. 

4.1 The Nuclear Modification Factor, Raa 

In Fig. 13, we plot Raa vs. px at energies (a) ^/snn = 22.5 GeV, (b) ^/snn = 62.4 GeV and 
(c) -\/snn = 200 GeV. The solid lines in the figures show the SCM-based results, wherein the 
experimental results are taken from Refs. [23], [25] and [26]. The dotted lines in the figures 
represent the pQCD-based calculations [17], [28]. Moreover, in Fig. 13(d) we have plotted the 
< Raa > vs. N par t between the range 2.5 < px < 3.5 Gev/c for the same collision and at the 
stated energies. Data-points are taken from the experiments by PHENIX Collaboration [17] . 
The solid lines in the Fig. 13(d) represent the our calculationl results and the dashed lines show 
the corresponding pQCD-oriented theoretical calculations made by Vitev [17] . [28] . 

4.2 The Nuclear Modification Factor, R C p 

In Fig. 14(a) we have plotted Rcp for the centralities — 10% and 60 — 94% against px for the 
7r°.The solid line in that Figure represents the SCM-predicted results arising out of the eqn.(25) 
and the dotted line in that gives the prediction from Quark Combination Model [27]. No data 
on Rcp for production of neutral pions at RHIC energies have yet been reported. 

Similarly, in Fig. 14 (b) we have plotted Rcp for the centralities — 10% and 40 — 60% 
against px for the K°. The solid line in the Figure 14(b) depicts the SCM-based plot against 
the experimental result [30] while the dotted line in the same figure represents the results attained 
by the Quark Combination Model [27] . 

5 Concluding Remarks 

The model applied here gives fair description of the px-spectra of all the light secondaries with 
the chosen values of the two parameters. Besides, the centrality-dependence of the pr-spectra 
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is also well-reproduced as is indicated by the figures. Slight disagreements observed at very 
Iow-pt (pt << 1 GeV/c) are due to the fact that the model has turned effectively into a mixed 
one with the entry of a power-law form due to the physics of partonic rearrangement factor. 
This dominance of power-law form disturbs, to a considerable extent, the agreement between 
data and model-based calculations for the extremely 'soft' (very low-py) values. Among the 
secondaries produced, the particle-antiparticle ratios and the proton-pion ratios are also in good 
agreement with the measured values. The obtained nuclear modification factors represented 
by Raa and Rcp (central-to-peripheral) are also in accord with the measurements. It is to be 
noted that we achieve all these with a new mechanism and introduction of some simple and basic 
ansatz like, the physics of large-pT nucleus-nucleus collisions and by introducing the properties 
of factorization, scale-breaking, mixed models with the combine of power-and-exponential laws, 
along with the principles of structural rearrangement factors at large transverse momenta. 

Selected comparisons of our model-based results with two other model-dependent calculations 
reveal neither sharp disagreement with any of them, nor very splendid agreement with either of 
them which are generically of standard model variety. Rather, on an overall basis, our results are 
in better agreement with data than either of them. This is a factor which is of some substance 
and importance to us. 

From a careful scrutiny of the fit-parameters we discover the following properties of them: 
(i) the structural rearrangement factor is clearly centrality dependent; it increases very slowly 
with gradual rise in the peripherality of the collisions, (ii) Secondly, the coefficients of the p\ 
in the exponential term are clearly energy-dependent in nature; on the contrary they manifest 
themselves to be independent of the centrality-measure of the interactions. 

The so-called suppression of the cross-sections at large-pr in heavy ion collisions is addressed 
here without resorting to the ideas of the 'jet-quenching' which is perceived to be to be one of 
the main conceptual pillars of heavy ion physics. 

An interesting question crops up in this connection. Here, we dealt with some aspects of 
Cu + Cu collisions at RHIC energies. We have chosen to remain silent about the physics of 
'quark-gluon plasma' (QGP) formation. The perturbative quantum chromodynamics (pQCD) 
predicted the formation of quark-gluon plasma (QGP). But the RHIC experiments failed to 
detect any plasma state; rather they found a "new kind of fluid state with very low viscosity" . 
So QCD prediction faltered at the first pillar. Secondly, how perfect the fluid observed at RHIC 
is cannot still be ascertained; the answer is not yet without caveats. Thirdly, we do not consider 
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this proposed 'plasma'-state to be any startling revelation, because when heated to very high 
temperatures, caused by the thermal motion of the molecules in the macroscopic matter, the 
solid substances melt down and turn into a variety of liquid. Almost in a similar manner, 
if microscopic particle-constituent matter is raised to very high temperatures attained by the 
extremely energetic collisions, the microscopic matter might also be converted to a liquid of 
somewhat unknown nature, and thus obviously to a 'new' kind |31j-|33j. So we do not pay much 
attention to the pQCD-based predictions on QGP and /or of suppression phenomena. 

Finally, we sum up by stating that (i) the model under consideration here explains and 
accommodates quite well the data on Cu + Cu collisions at various energies; and (ii) Quite 
agreeably, the values of Rcp obtained by the present calculations are not in good agreement 
with data. This could be attributed to our neglect of the effects of final state re-scatterings and 
some other complex physical factors. 
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Table 2: Values of a, Nr and /3 for tt and tt + productions in Cu + Cu collisions at v /J/v r /v=22.5, 
62.4 and 200 GeV 



Centrality 



fs]^=22.b GeV 





TT 

an- Nf /V 


7T+ 

a n + Nf 


0-10% 
10-30% 
30-60% 
60-100% 
Minbias 


0.901 3.454 0.703 
0.492 3.444 0.703 
0.239 3.413 0.703 
0.033 3.410 0.703 
0.395 3.431 0.703 


1.001 3.454 0.703 
0.572 3.444 0.703 
0.244 3.413 0.703 
0.034 3.410 0.703 
0.468 3.431 0.703 


Centrality 


TT~ 

a n - Nf P n - 


.2.4 GeV 

7T+ 

OV+ Nf f3 n+ 


0-10% 
10-30% 
30-60% 
60-100% 
Minbias 


0.707 3.035 0.468 
0.464 3.030 0.468 
0.180 3.026 0.468 
0.038 3.016 0.468 
0.204 3.142 0.468 


0.786 3.035 0.468 
0.474 3.030 0.468 
0.228 3.026 0.468 
0.044 3.016 0.468 
0.207 3.142 0.468 







Nf 


TT 

Pn- 


7T+ 

a n + 


Nf 




0-5% 


1.098 


3.597 


0.293 


1.113 


3.597 


0.293 


5-10% 


0.898 


3.572 


0.293 


0.910 


3.572 


0.293 


10-15% 


0.750 


3.567 


0.293 


0.760 


3.567 


0.293 


15-20% 


0.687 


3.552 


0.293 


0.696 


3.552 


0.293 


20-30% 


0.635 


3.544 


0.293 


0.643 


3.544 


0.293 


30-40% 


0.523 


3.525 


0.293 


0.530 


3.525 


0.293 


40-50% 


0.353 


3.518 


0.293 


0.358 


3.518 


0.293 


50-60% 


0.283 


3.508 


0.293 


0.287 


3.508 


0.293 


60-70% 


0.113 


3.491 


0.293 


0.115 


3.491 


0.293 


70-80% 


0.085 


3.476 


0.293 


0.086 


3.476 


0.293 


80-92% 


0.004 


3.453 


0.293 


0.004 


3.453 


0.293 
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Table 3: Values of a, Nr and /3 for K and K + productions in Cu+Cu collisions at ^/snn=22.5, 
62.4 and 200 GeV 



Centrality 






K~ 


52.5 GeV 
K+ 








a K - 




Pk- 


a K + 






0-10% 


0.406 


2.304 


0.863 


0.469 


2.304 


0.863 


10-30% 


0.185 


2.285 


0.863 


0.255 


2.285 


0.863 


30-60% 


0.069 


2.274 


0.863 


0.114 


2.274 


0.863 


60-100% 


0.012 


2.270 


0.863 


0.018 


2.270 


0.863 


Minbias 


0.118 


2.264 


0.863 


0.151 


2.264 


0.863 


Centrality 






y/SNN=t 

K~ 


S2.4 GeV 
K+ 








a K - 




Pk- 


a K + 


Nf 




0-10% 


0.675 


2.714 


0.571 


0.755 


2.714 


0.571 


10-30% 


0.405 


2.704 


0.571 


0.524 


2.704 


0.571 


30-60% 


0.162 


2.688 


0.571 


0.187 


2.688 


0.571 


60-100% 


0.023 


2.658 


0.571 


0.028 


2.658 


0.571 


Minbias 


0.217 


2.681 


0.571 


0.251 


2.681 


0.571 


Centrality 






\/snn=- 
K~ 


200 GeV 

K+ 








a K - 




Pk- 


a K + 


Nf 


Pk+ 


0-5% 


0.286 


2.939 


0.417 


0.294 


2.939 


0.417 


5-10% 


0.250 


2.839 


0.417 


0.257 


2.839 


0.417 


10-15% 


0.214 


2.819 


0.417 


0.220 


2.819 


0.417 


15-20% 


0.210 


2.805 


0.417 


0.216 


2.805 


0.417 


20-30% 


0.181 


2.795 


0.417 


0.186 


2.795 


0.417 


30-40% 


0.165 


2.765 


0.417 


0.170 


2.765 


0.417 


40-50% 


0.112 


2.734 


0.417 


0.115 


2.734 


0.417 


50-60% 


0.051 


2.714 


0.417 


0.053 


2.714 


0.417 


60-70% 


0.041 


2.698 


0.417 


0.042 


2.698 


0.417 


70-80% 


0.024 


2.672 


0.417 


0.026 


2.672 


0.417 


80-92% 


0.013 


2.652 


0.417 


0.014 


2.652 


0.417 
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Table 4: Values of a, Nr and (3 for p and p productions in Cu + Cu collisions at v /J/v r /v=22.5, 
62.4 and 200 GeV 



Centrality 






^s NN =22.b GeV 










p 

ap (3 p 


V 

a p N P R f3 p 




0-10% 


0.211 


0.826 


0.853 


0.344 


0.826 


0.853 


10-30% 


0.111 


0.780 


0.853 


0.129 


0.780 


0.853 


30-60% 


0.044 


0.695 


0.853 


0.051 


0.695 


0.853 


60-100% 


0.008 


0.691 


0.853 


0.010 


0.691 


0.853 


Minbias 


0.070 


0.798 


0.853 


0.101 


0.798 


0.853 


Centrality 






^s NN =V2A GeV 










p 

a p Nft (3 P 


V 

a p N P R f3 p 




0-10% 


0.468 


1.118 


0.618 


0.957 


1.118 


0.618 


10-30% 


0.256 


0.958 


0.618 


0.513 


0.958 


0.618 


30-60% 


0.108 


0.945 


0.618 


0.192 


0.945 


0.618 


60-100% 


0.012 


0.930 


0.618 


0.023 


0.930 


0.618 


Minbias 


0.184 


1.178 


0.618 


0.298 


1.178 


0.618 


Centrality 






^/s NN =20U GeV 










p 

a p N R (3 p 


V 

a p (3 P 




0-5% 


0.166 


1.251 


0.426 


0.185 


1.251 


0.426 


5-10% 


0.156 


1.231 


0.426 


0.181 


1.231 


0.426 


10-15% 


0.113 


1.211 


0.426 


0.131 


1.211 


0.426 


15-20% 


0.105 


1.192 


0.426 


0.122 


1.192 


0.426 


20-30% 


0.087 


1.172 


0.426 


0.101 


1.172 


0.426 


30-40% 


0.074 


1.154 


0.426 


0.086 


1.154 


0.426 


40-50% 


0.038 


1.133 


0.426 


0.044 


1.133 


0.426 


50-60% 


0.028 


1.112 


0.426 


0.032 


1.112 


0.426 


60-70% 


0.017 


1.105 


0.426 


0.020 


1.105 


0.426 


70-80% 


0.006 


1.099 


0.426 


0.007 


1.099 


0.426 


80-92% 


0.003 


1.075 


0.426 


0.004 


1.075 


0.426 
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p T (GeV/c) p T (GeV/c) 

(a) (b) 



p+p =k* +X (s NN 1/2 =200 GeV) 




1e-005 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

p T [GeV/c] 



(c) 

Figure 1: Plots for ir production in p + p collisions at energies (a) ^/snn = 20 GeV, (b) v / s]v]v = 63 GeV 
and (c) y^SNN = 200 GeV. Data are taken (a) from Ref. [12], (b) from Ref. Q]5] and (c) from Ref.s [T3j . 
[IS] , Solid lines in the Figures show the SCM-based plots. 
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Cu+Cu=7i"+X[(s NN ) =22.5 GeV] 
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Figure 2: Centrality dependence of the pt distribution for ir for different centralities and at energies 
(a) 22.5 GeV [TJ], (b) 62.4 GeV 18, and (c) 200 GeV 19, in Cu + Cu collisions. The solid lines in the 
Figures 2(a), 2(b) and 2(c) show the SCM calculations for different centralities. 
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Figure 6: Centrality dependence of the px distribution for p for different centralities in Cu+Cu collisions 
at (a) ^/IW^ =22.5 GeV [TS], (b)y/s NN =62.4 GeV Q2] and for (c) y/s NN =200 GeV IS]. The solid 
lines in the Figures show the SCM calculations . 
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Figure 7: Centrality dependence of the px distribution for p for different centralities in Cu+Cu collisions 
at (a) =22.5 GeV [TS], (b)y/s NN =62.4 GeV Q2] and for (c) y/s NN =200 GeV IS]. The solid 

lines in the figures show the SCM calculations . 
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Figure 8: The N part versus it /ir + ratio behaviours at (a) ^/sjvjv =62.4 GeV and (b)y/s NN =200 GeV. 
Data are taken from PHOBOS [21], [22] ■ The solid line in Fig. 8(a) and the filled circles in Fig. 8(b) 
show the theoretically calculated values. 
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Fi gure 9: The Np ar t versus K /K^ ratio behaviours at (a) -\/snn — 62.4 GeV and (b)-\Zsjyjy =200 GeV. 
Data are taken from PHOBOS [21], [22]. The solid line in Fig. 9(a) and the filled squares in Fig. 9(b) 
show the SCM-based calculated values. 
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Figure 10: The p/p ratio behaviours as a function of number of participants (N part ) in Cu + Cu reaction 
at (a) ^/IWn =62.4 GeV and (b) at ^i^F =200 GeV. Data are taken from PHOBOS [21], [22]. The 
solid line in Fig. 10(a) and the filled circles in Fig. 10(b) show the SCM-based calculated values. 




Figure 11: Ratios of (a) p/tt" and (b) p/ir + as a function pr for central Cu + Cu reactions at ^/snn 
=22.5, 62.4 and 200 GeV. Data in these Figures are taken from [23] , The solid lines show the SCM-based 
results. 
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Figure 12: Comparisons (a) between the SCM, the Quark Combination Model [27] and the pQCD 
approach [55] for the tt° production at 200 GeV[T2] and (b)between the SCM, the Quark Combination 
Model [27] for the Kg production for Cu + Cu reactions at ^snn =200 GeV. Data are taken from [29] 
and [3D]- The solid lines in those figures represent SCM-based results, wherein the dashed and dotted 
lines represent the Quark Combination Model [27] and the pQCD approach [28 . 
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Figure 13: Plots of px versus Raa (as defined in the text) at energies (a) ^snn =22.5 GeV, (b) ^snn 
= 62 .4 and (c) ^snn = 200 GeV. Data in these Figures are from [23], [25] and [26]. The solid lines show 
the SCM-based results, wherein the dashed lines represent the pQCD-oriented calculations [T7],[28]. (d) 
Plot of < Raa > vs. N part for pr-ranges like 2.5 < pr < 3.5 GeV/c. Comparisons of the nature of 
average nuclear modification factors based on two sets of calculations, one done by SCM and represented 
by solid lines and the other (dashed lines) obtained and shown by Vitev [17], [28] are made here. 
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Figure 14: Plots for the Rcp behaviours versus pt of (a) tt° and (b)K°. Pion data for Cu + Cu collisions 
has not yet available. Plots in Fig. 14(a) are predictive comparison between SCM-based result and result 
from Yun-Fei et al.[27]. Data in Fig. (b) are taken from STAR [3U]- The solid line in Fig. 14(b) shows 
the SCM-based results while the dotted line represents the results of Yun-Fei et al.[27]. 
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